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Abstract. Dynamical systems methods are used to investigate global behaviour of 
the spatially flat Friedmann-Robertson-Walker cosmological model in gravitational the¬ 
ory with a non-minimally coupled scalar field and a constant potential function. We 
show that the system can be reduced to an autonomous three-dimensional dynamical 
system and additionally is equipped with an invariant manifold corresponding to an 
accelerated expansion of the universe. Using this invariant manifold we find an ex¬ 
act solution of the reduced dynamics. We investigate all solutions for all admissible 
initial conditions using theory of dynamical systems to obtain a classification of all 
evolutional paths. The right-hand sides of the dynamical system depend crucially on 
the value of the non-minimal coupling constant therefore we study bifurcation values 
of this parameter under which the structure of the phase space changes qualitatively. 
We found a special bifurcation value of the non-minimal coupling constant which is 
distinguished by dynamics of the model and may suggest some additional symmetry 
in matter sector of the theory. 
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1 Introduction 

In the modern cosmology there is an increasing interest in the description of a matter 
content of the universe in terms of a scalar held. The accelerated expansion of the 
universe [1, 2] stimulated investigations in the held of dynamical dark energy models 
[3]. The quintessence idea was formulated [4, 5] as a simplest model involving a scalar 
held with a potential function in order to describe the current accelerated expansion of 
the universe. The simplest candidate for the dark energy driving the current phase of 
the evolution of the universe seems to be a positive cosmological constant, however such 
a explanation suffers from the fine tuning problem [6, 7] and the coincidence problem 
[7, 8]. Nevertheless the straightforward application of this assumption within the 
general theory of relativity together with the Friedmann-Robertson-Walker symmetry 
leads to the ACDM model which is a standard model of the cosmological evolution 
[9, 10] and is favoured by the observational data [11-13]. Number of alternatives were 
proposed in order to alleviate and eliminate problems with a cosmological constant 
term like a phantom dark energy [14, 15] or a extended quintessence [16-22], 

The inclusion into a matter sector of the gravitation theory of a non-minimal 
coupling term between the gravity and the scalar held of the form — £R(p 2 [23-25] leads 
to natural and simplest generalisation of the scalar held Lagrangian. The motivations 
for this term come form different contexts. The general relativity has the method¬ 
ological status of an effective theory with a given Lagrangian and such contribution 
naturally emerge in its expansion [26]. The non-minimal coupling between the scalar 
curvature and the scalar held appears as a result of quantum corrections to the scalar 
held in curved space [27-29] and is required by the renormalisation procedure [24], The 
non-minimal coupling is also interesting in the context of superstring theory [30] and 
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induced gravity [31]. While the simplest inflationary model with a minimally coupled 
scalar field and a quadratic potential function is no longer favoured by the observational 
data [32-34] there is an urgent need to extend this paradigm further beyond the mini¬ 
mally coupled case. From the theoretical point of view the coupling constant becomes 
a free parameter in the model which value should be estimated from the observational 
data [35-38] or should be obtained from a more fundamental theory. The model with 
a zero value of the non-minimal coupling was investigated in detail by Faraoni [39]. In 
this case Faraoni have found an exact solution which can be important in the context 
of an inflation description. Our investigation generalises Faraoni’s results to the case 
of the constant non-minimal coupling. 

The non-minimally coupled scalar field cosmology was investigated by many au¬ 
thors in the connection with an inflationary epoch as well as a description of the 
quintessence [20, 40-63]. Uzan [52] pointed out the significance of scaling solutions in 
generalised inflation scenarios. The question of integrability of non-minimally coupled 
scalar field cosmological models has been recently studied by Kamenshchik et al. [64] 
who looked for integrable cases with a different forms of scalar field potential function. 
The new dynamical system methods has been recently developed by Skugoreva et al. 
[65] in the context of a more general form of coupling parametrised by the scalar field. 
In the standard model of particle physics a non-minimally coupled Higgs field plays 
also important role [66-69]. 

The dynamical systems methods are widely used in cosmological applications since 
seminal papers by Belinskii [70, 71] and the most widespread parameterisation of the 
phase space is by the so-called expansion normalised variables [72] (see collection of 
works on the dynamical system analysis of anisotropic models [73]). 

In the present paper we investigate a dynamical behaviour of the flat Friedmann- 
Robertson-Walker cosmological model with matter in the form of the non-minimally 
coupled scalar field with a constant potential function and a barotropic dust matter. 
This model is the simplest generalisation of the ACDM model and can be treated as 
a testing ground for evolutional paths of the universe. Adding the effects of the non- 
minimal coupling to the general relativity models exact solutions can be found only 
for special cases of the non-minimal coupling, namely for the minimal £ = 0 and the 
conformal £ = |. The dynamical system is integrable and the exact solution can be 
obtained in the form of elliptic functions like for the standard Friedmann-Robertson- 
Walker models with dust, cosmological constant and radiation matter, in the case of 
the conformal coupling, or a stiff matter contribution in the case of minimal coupling. 
While there is no possibility to obtain exact solutions in general case one can apply 
the dynamical system methods for studying all evolutional paths, i.e. trajectories 
in the phase space. We demonstrate that for the flat Friedmann-Robertson-Walker 
model the dynamics can be reduced to the form of the three dimensional autonomous 
dynamical system in which the coupling constant becomes the free parameter. The 
main advantage of the dynamical system methods is that that without knowledge of 
an exact solution it is possible to investigated the properties of the solutions as well as 
their stability. If the right sides of the system vanish then we obtain solutions in the 
form of critical points representing, from the physical point of view, asymptotic states 
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of the system. In the phase space, which is the space of states of the system, they are 
represented by critical points of different types. Global dynamics is represented by the 
global phase portrait which contains all information about behaviour of trajectories 
and stability of the asymptotic states. Important information about the behaviour of 
the system is encoded in the behaviour of trajectories on invariant manifolds of lower 
dimension, which we identify for the system and we find exact solutions for the reduced 
dynamics. 

The complete analysis of the dynamics requires the investigation of trajectories at 
the infinity. For this aim the construction of the Poincare sphere with a corresponding 
projection can be useful (see [74, 75]). An alternative construction which can be applied 
in the study of dynamics is to introduce the projective coordinates on the phase space. 
Then projective maps cover the dynamics at infinity. Such an analysis of the behaviour 
of trajectories at infinity is important if we are going to formulate of some conclusions 
about the global stability of the system. 

Our investigation explicitly shows a degree of complexity of the dynamics as 
the value of the non-minimal coupling constant changes. We construct the phase 
portraits for all representative cases of the non-minimal coupling constant. The global 
dynamics crucially changes if this constant assumes some critical value - a bifurcation 
value. Therefore the global phase portraits offer us classification of a scalar field 
cosmology with the constant value of the non-minimal coupling and can be treated as 
generalisation of corresponding classification of the ACDM model. Of course due to 
conformal duality of the system under considerations to the minimally coupled case, 
any change of the non-minimal coupling parameter £ will lead to the change of the 
corresponding non-constant potential in the Einstein frame formulation of the theory. 
One can admit that more complicated potential functions lead to more complicated 
dynamics. In the Jordan frame, followed in this paper, we End that the degree of 
complexity changes with the value of the non-minimal coupling constant with respect 
to a simple minimally or conformally coupled case which from the other hand lead to 
simple integrable dynamics. 

The analysed dynamical system poses a three dimensional phase space. In the case 
of a zero potential (the case without the cosmological constant) the right hand sides of 
the system are homogeneous functions of degree three which unable us to reduce the 
system to the two-dimensional one. In the three-dimensional phase space the dynamics 
is very complex and crucially depends on the constant value of the coupling. We point 
out that the value of £ = 3/16 is distinguished by the dynamical behaviour of the 
model and leads to physically interesting and appealing evolution of the universe. 

2 The model 

The working model consists of the spatially flat Friedmann-Robertson-Walker (FRW) 
universe filled with the non-minimally coupled scalar field and a barotropic matter 
with the constant equation of state parameter w m . The total action is in the following 
form 

S = S g + S <p + S m , (2.1) 
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where the gravitational part of the theory is described by the standard Einstein-Hilbert 
action integral 

Ss = hI d ‘ x ^ R ’ < 2 - 2 > 

where k 2 = 87 tG, and the matter part of the theory is described by the following action 
for a non-minimally coupled scalar held 

S^-IJ d [eV a <f> V Q 0 + e^Rcj) 2 + 2f/(0)) , (2.3) 

where e = + 1,-1 corresponds to the canonical and the phantom scalar held, respec¬ 
tively, and a barotropic matter 


S m = J d 4 x^gC m . (2.4) 

The held equations for the theory are 

Rf* - \s^R = K 2 (TjS + T«) , (2.5) 

where the energy-momentum tensor for the non-minimally coupled scalar held is given 
by 


TjS = - e\g^V a <j>V a <i> - f/(0)^+ 

l ( 2 - 6 ) 

+ {R^lu - 2 9tiuR)4> 2 + j , 

and for the barotropic matter included in the model 

^ (Pm T v T PmdfJ.i' • (2-7) 

Let us note that when the non-minimal coupling between the scalar held and 
the curvature is present the held equations can be written in several nonequivalent 
ways. In the case adopted here the energy momentum tensor for the scalar held 
(2.6) is covariantly conserved, which may not be true for other choices [76-78]. For 
example, the redefinition of gravitational constant = k~ 2 — e £</> 2 makes it time 
dependent. The effective gravitational constant can diverge for the critical value of the 
scalar held (f> c = ±(e£ft 2 )~ 1//2 . Though the FRW model remains regular at this point, 
the model is unstable with respect to arbitrary small anisotropic and inhomogeneous 
perturbations which become infinite there. This results in the formation of a strong 
curvature singularity prohibiting a transition to the region k 2 s < 0 [79, 80]. 

In cosmological applications one can obtain substantial mathematical simplifica¬ 
tion using conformal transformation techniques, especially in the absence of ordinary 
matter, in order to relate cosmological models with a non-minimally coupled scalar 
held with its conformal counterpart with a minimally coupled held. The Jordan frame 
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action integral, where the scalar held is non-minimally coupled to the Ricci scalar cur¬ 
vature is mapped into an Einstein frame where now the transformed scalar held 0 is 
minimally coupled. The two frames are physically nonequivalent unless variable units 
of time, length, and mass are adopted in the Einstein frame [77, 78]. 

In the theory under considerations using the conformal factor fl 2 = 1 — £^K 2 (f) 2 
the Einstein frame counterpart of the action integrals (2.2) and (2.3) is 


q(E) _ 
i3 9+4> 


1 

2 


d 4 x 


v^l 


( J_6 _ 1 -g^(! - 6Q/cV 

K 2 £ (1 - £^K 2 (j) 2 ) 2 




n UM \ 

(1 - £^ K 2 ( j ) 2 ) 2 ) ’ 


and the scalar held can be transformed as 


- \A(! -efC 1 - 6 O« 2 0 2 ) 

# = 1 -7-^' 

1 — £^K z (p z 

The condition for the theory given by this action integral to be ghost free is 

e(l -ef(l - 6O« 2 0 2 ) > 0, 

provided that fl 2 = 1 — £^n 2 (j) 2 > 0, in this sense the conformal transformation tech¬ 
nique cannot provide solutions with the scalar (j) crossing the barriers (f> c = ±(e^K 2 )^ 1 / 2 
and has the same limitations of the form of the held equations using the effective cou¬ 
pling = k~ 2 - £^(p 2 [77, 78]. 

In the present paper we work exclusively in the Jordan frame formulation of 
the theory leaving aside the contentious question concerning the physical equivalence 
between the Jordan frame and the Einstein frame [81, 82], 

From now on we assume that the barotropic matter in the form of a dust matter 
and taking the spatially hat FRW metric 

ds 2 = —dt 2 + a 2 (t ) ^dx 2 + d y 2 + d z 2 ^ , 

we obtain the energy conservation condition in the standard form 

~~^H 2 — Pc/, + p m (2-8) 

where p m is the energy density of the barotropic dust matter and 

Pc/, — £-</> 2 + Uo + £3^H 2 (j) 2 + £6 f^Hcfxj) 

is the energy density of the non-minimally coupled scalar held and we have assumed a 
constant potential function U(4>) = Uq = const. The dynamical equation for the scalar 
held is 

4> + 377 (j) + £R(j) = 0 . 

Now, we divide the constraint equation (2.8) by the present value of the Hubble 
function Hq = H(ao) and we obtain 

^ + fl TO ,o + e(l — 6£)x 2 + e6£(x + z) 2 (2.9) 
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where 


= k 2 U 0 _ n 2 p mfi 

- 3H 2 , m ,0 - 3H 2 

and we introduced new dimensionless variables, which are the energy phase space 
variables in a slightly different form, 


x = 


k c f> 


z = 


K H 
Vd Hq 


( 2 . 10 ) 


The acceleration equation expressed in the new variables is 

H ^(1 — 6£)x 2 + 212 a,0 + 

H2 ~ + 

and using the energy conservation condition one can eliminate the h2 mi o term. 

The dynamical system on variables x, z and h — M- is in the following form 


da; 

din a 


—3a; — 6 £z 




dz H 

— = X + Z 1H' 
dm a H 2 

dh _ H 

dh^ _ 'IP ’ 


( 2 . 11 ) 


where 

H 0 3 ^a,o + \h 2 — e(l — 6^)x 2 — e2£(x + z) 2 

H 2 2 h 2 - e6^(l - 6 £)z 2 ' 1 ' ' 

From the last equation in the dynamical system (2.11), one can easily notice, that 
the system has two invariant manifolds h = H / Hq = 0 and H / H 2 = 0. This second 
property will be used later in the paper in order to obtain an explicit solution of a 
restricted dynamical system. 

We clearly see that the right hand sides of the dynamical system (2.11) are rational 
functions of their arguments due to form of the acceleration equation (2.12). To remove 
this degeneracy, we perform the following time transformation 

and the dynamical system describing the model under investigations is now in the 
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following form 


—- = -3 x(h 2 - e 6 £(l - 6 £) 2 2 ) — 

- 9fz(fi A) o + \h 2 - e(l - 60a; 2 - e2£(x + z) 2 ^ , 


—- = x(h 2 — e6£(l — 6£)£ 2 ) + 
dr v ' 

+ ^2 ^a,o - h 2 - e(l - 60a; 2 - e2£(x + z) 2 + ^8£(1 - 6£)z 2 ^ , 


(2.14) 



where all functions on the right hand side are polynomials and the time transformation 
(2.13) removes singularities of the system (2.11), where h 2 — e6£(l — 6£)z 2 = 0. 

The dynamical system (2.14) holds a very special property. For the vanishing po¬ 
tential of the scalar held f^o = 0 the right hand sides of this system are homogeneous 
polynomial of degree 3 in the phase space variables. In such case the projective coor¬ 
dinates can be used in order to lower dimension of the dynamical system and decouple 
one of the equations [74, 75]. This property will be widely used in subsection 3.2 while 
analysing dynamical behaviour of the system at infinity of the phase space. 

Before we proceed to dynamical analysis of the model, we need to consider two 
special values of the non-minimal coupling constant £, namely, minimal £ = 0 and 
conformal coupling £ = | values. 

Let consider first the case of the minimal coupling. For £ = 0 from (2.9) we have 



(2.15) 


and the system (2.11) reduces to single equation 



which can be easily integrated resulting in 



where the initial condition x(ao) is taken at the present epoch. Now, inserting this 
solution in to (2.15) we obtain the following Hubble function 



(2.16) 


which describes the background dynamics of the spatially flat FRW model with the 
minimally coupled scalar held. The last term indicated contribution from the scalar 
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field and has a similar form as a stiff matter with the equation of state parameter 
w m = 1. Depending on the type of scalar held (canonical e = +1 or phantom e = —1) 
this contribution has positive or negative energy density. The observational constraints 
on the such modified Hubble function were previously investigated in [83] leading to 
the conclusion that this model is indistinguishable from the ACDM model. 

The second case considered is the conformal coupling. For £ = | from (2.9) we 
have „ 

— ^a,o + ^m,o (—+ e (x(a) + z(a)) 2 (2-17) 

£=i \ a o/ 

and the system (2.11) reduces to 



d(x + z) 
din a 


—2(x + z ), 


which can be explicitly integrated resulting in 


x(a) + z(a) 


(x(a 0 ) + z(a 0 )) 
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where the initial conditions x(ao) and z(ao) are taken at the present epoch. Inserting 
this solution in to (2.17) we obtain 


H(a ) 

H(ao) 


2 



f^A,o T o 



-3 

+ £ {xo + zo ) 2 



(2.18) 


The resulting Hubble function describes the standard ACDM model with additional 
term coming from the conformally coupled scalar held which behaves like a radiation 
in the model. The energy density of this radiation-like term is positive for a canonical 
scalar held and is negative for a phantom scalar held. Recently, Boisseau et al. found 
a bouncing solution with conformal coupling tending to dS space and containing dark 
radiation [84], 

Using specihc assumptions regarding vanishing of the scalar held potential func¬ 
tion or the barotropic matter in the model the both of Hubble functions (2.16) and 
(2.18) can be integrated in exact form in terms of elliptic functions [85, 86]. 


3 The non-minimal coupling and the structure of the phase 
space 

3.1 Critical points and stability 

In dynamical system theory phase space of a system is organised by asymptotic states, 
i.e. critical points and trajectories joining them. In section 2 we were able to reduce 
dynamics of the model under investigations to the 3—dimensional dynamical system 
(2.14) where right hand sides are polynomials in the phase space variables. Redefini¬ 
tion of a time variable made possible to move from rational functions to polynomial 







functions of the phase space variables and remove singularities of the dynamics. The 
price to pay is a new time function along the phase space trajectories. 

In this subsection we are interested in a local behaviour of the dynamics which we 
want to directly connect with the Hubble function describing the background evolution 
of the cosmological model. The most important structures in the dynamical system 
theory are asymptotic states of the dynamics in the form of critical points. Our analysis 
we begin from the system (2.11) with the natural logarithm of the scale factor as time 
variable along the phase space trajectories 


dx 

din a 


—3a; — 6 




dz H 

dm a H z 

dh _ H 
dh^ “ l ~FP ’ 


(3.1) 


where 

H _ dVt Kfl + \h 2 -£(l-d£)x 2 -£2 £(x + z) 2 

H 2 ~ + 2 h 2 - (1 - 6 Oz 2 ' 1 j 

In what follows we assume that the non-miuimal coupling constant £ is different 
from the minimal and conformal coupling values £ ^ 0, | and that the phase space 
variables £ and h are such that there is no singularity in the denominator of the 
acceleration equation (3.2). 

Dynamical system (3.1) admits four non-singular critical points at a finite domain 
of the phase space. Two of them located at 

{x* = 0,z* = 0,h* = ±y/n^o), (3.3) 


correspond to the de Sitter and anti-de Sitter states respectively, with the vanishing 
acceleration equation (3.2). The only condition for the existence of those points is 
Da,o > 0 and it does not depend on a value of the non-minimal coupling constant £. 
The remaining two asymptotic states located at 


x 


= T4£i 


n 


A,0 


£ 2 ^( 1 - 60 ( 3 - 100 ’ 


= ±( 1 - 20 - 


n 


A,0 


£ 201 - 60 ( 3-100 


,h* = 0 


(3.4) 

correspond to the Einstein static solutions with the vanishing Hubble function and 
the asymptotic value of the effective equation of state parameter w e g = —1 — 

The conditions for existence are — , 2 qi-6‘^'':-ioo > ^ together with the non-minimal 
coupling constant different form £ 0 {0, T}. Explicitly this leads to the following 
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conditions for values of the cosmological constant and the non-minimal coupling 


^A,o > 0 : 

£ 

= -l, 

0 < £ < 1 

V 


£ > — 

10 ’ 


£ 

— +1 j 

£ < 0 

V 

1 

6 

A 

<r r> 

A 

ol“ 

^A,o < 0 : 

£ 

= -l, 

£<o 

V 

i 

6 

< £ < — 
^ ^ ^ 10 


£ 

— +1 7 

0<£<i 

V 


£ > — 

^ ^ 10 • 


Now we can proceed to finding linearised solutions in the vicinity of those points 
and to obtain a physical interpretation of the dynamics. The de Sitter and anti-de 
Sitter states are given by the following critical points 

(x* = 0,zr = 0,(h*) 2 = Q Afi ), (3.6) 

in terms of the original scalar held variables it is 

(£>* = 0, 4>* = 0 , H* = const.) . 

The linearisation matrix calculated at this point lead to the following eigenvalues 

Aig = — 2 T 2 \J^- —> A 3 = -3. (3.7) 

The stability condition for this critical point is that the real parts of Ai and A 2 eigenval¬ 
ues must be negative Re(Ai) < 0 and Re(A 2 ) < 0. This leads to the following condition 
£ > 0, the non-minimal coupling constant £ must be positive in order to obtain the 
asymptotically stable de Sitter state. In details we have 


a stable node : 0 < £ < ^ , 

a stable focus (a sink) : £ > , 


(3.8) 


while for £ < 0 the de Sitter state is a transient state of a saddle type critical point 
with Ai < 0, A 2 > 0 and A 3 < 0. The linearised solutions to the dynamics in the 
vicinity of the de Sitter and anti-de Sitter states are 


x(a) = ^hr 2 ( (Ai Ax ~ l2 ^ z) ~ (A2 Al “ 12? Az) ( i»T) • 

z(o) = yvy ( (Aa: - A2Az) (Jr)"' - (Al - Ai Az) {£<>T ) ■ 

h(a) = h," + Ah (J>) . 


(3.9) 


where Ai and A 2 are the two out of three eigenvalues of the linearisation matrix and 
Ax = xS l \ Az = z® and Ah = h® — h* are the initial conditions. 
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From the last equation of this system we directly obtain that up to linear terms 
in initial conditions the Hubble function is 

2 / \ -3 


(H(a)\ 

W<*> )J 


^A,0 “1“ ^m,0 


(3.10) 


and is indistinguishable from the Hubble function for the standard ACDM model. 
The critical points corresponding to the static Einstein solutions are give by 


x = 


4£ 


2 *, oo = -e(i - 2 zy 


^A,0 


l-2r ,vw 2^(1-60(3-100 

and in the original scalar held variables this is equivalent to 

(</>* = const., 4>* = ±oo , H* = 0). 

The linearisation matrix gives the eigenvalues 

2 . 4£ 


,h* = 0 


(3.11) 


Ai — —3 , A 2 — —5 + 


At — 


(3.12) 


1 - 2 £ ’ 1-2C 

which constitute stability conditions for the state. Matching the existence conditions 
(3.5) we obtain the following stability conditions for a stable node : Ai < 0, A 2 < 
0 , A 3 < 0 


(3.13) 


^A,o > 0 : 

e = -1, 



£ = +1 , 

e<0, 

A1a,o < 0 : 

£ = - 1 , 

e<o, 


£ = +1 , 

e>|, 

for a saddle : Ai < 0 , A 2 > 0 , A 3 > 0 



A1a,o > 0 : 

e = -!, 

— < £ < 1 

^A,o < 0 : 

£ — +1 , 

A < £ < 1 
10 ^ ^ ^ 2 

and for a different type of saddle : Xi 

< 0 , A2 < 0 , A3 > 0 

^A,o > 0 : 

e = -!, 

°<£<! 


£ — +1, 

1 < f < A 
6 ^ 10 

^A,o < 0 : 

£ = ~l, 

1 <y < A 
6 ^ ^ ^ 10 


£ = +1, 

o<e<| 


Finally, the linearised solutions in the vicinity of the critical point corresponding to 
the Einstein static solution are 

1 


x(a) =x*~ ^(Ax + 3Az) + |(Ax + Az) , 


a \ x 2 


z(a = z 


+ J(Ax + 3Az) (Jy) 1 - 7 ;(Ax + Az) 


a \ x 2 


VaW 


(3.14) 


h (a) = &h (Ay' ’ 
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where x* and z* are the coordinates of the critical point, Ai, A 2 , A 3 are the eigenvalues 
of the linearisation matrix, and Arc = x^' — x*, Az = z^' — z*, Ah = h W are the initial 
conditions taken in the close vicinity of the critical point. 

From the last equation we directly obtain the evolution of the Hubble function 

H(a) _ H{a «) / a 
H(ao) ~ H(ao) 


which vanishes up to quadratic term in initial condition 


H(a) 

H(a 0 ) 


2 


0 . 


(3.15) 


The critical points investigated in this section constitute non-singular asymptotic 
states in the form of de Sitter state and static Einstein solution. Examination of the 
stability conditions (3.8) and (3.13) indicates that the both states can be stable during 
expansion only for the phantom scalar held £ — —l with f^o > 0 and £ > A. Most 
interestingly for the canonical scalar held e — +1 with Ha,o > 0 and negative values of 
the non-minimal coupling constant £ < 0 the de Sitter state is transient, in the form 
of a saddle type critical point, while the static Einstein state is asymptotically stable. 


3.2 The analysis of dynamics at infinity 

In order to obtain complete information about the structure of the phase space of the 
dynamical system (2.11) it is necessary to investigate the behaviour of the trajectories 
of this system at infinity. 

To study the dynamical system (2.11) we need to introduce three charts f7( n ) of 
projective coordinates 

l z h 

U W : U W = - ) V (1) = -» w W = - > 

U {2) : U( 2) = | , V(2) = \ , W(2) = ^ , (3-16) 

X z 

U {3) : u (3) = J x > m = ^ > w (3) = h ■ 

These three charts enable us to find a complete set of the critical points of the system 
(2.14) at infinity of the phase space. Some critical points will be present in more than 
one chart. 
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* 

U (D 

W (l) 

existence 

A 

1-4? 

2? 

0 

V£ e R \ {0} 

B 


0 

£ < 0 or £ > ± 

C 

1 


e = —1: £<0V£>|ore = +1: 0 < £ < \ 

D 

l± v /-2(l-6?) 

±\O 6 0! - 60^(1) 

e = - 1: £ > \ 

6? 

E 

0 

0 

V(eR 


Table 1. Critical points of the system (3.17) in the map £/)i) at infinity of the phase space 
defined by the condition = 0. 


The use of Um chart transforms the system (2.14) in to the following one 
^ = 3 u {1) ((' w 2 {1) - e6£(l - 60 ^( 2 i )) + 

+ 3£w(i) (^a,o u Ii) + — — — £ 2£(1 + W(i)) 2 j j , 

^ = (l + v ( i))(t«(i) -e6e(l-60v(i)) + 

+ t^OO^ 6 ^(i)) + -W(i) — e(l — 60 — e2^(l + wp)) 2 ^ , 

^ = w a )(Wd - e6 £(! - 6 04))+ 

+ “*■ 6^(1)) ^a,o m (i) + 2 ^( 1 ) — e (l — 60 — e2£(l + V(i)) 2 j j , 

(3.17) 


where the time function along the phase space trajectories is 


c%) 


= U, 


(1)^7 = 6 0«(i)) 


din a 


(3.18) 


The critical points of this system at the invariant manifold at infinity defined by the 
condition up) = 0 are gathered in table 1. 

The energy conservation condition (2.8) expressed in the projective coordinates 
of the chart C/p) is given by 


o T - -m.O 



M (1) = w (i) - £ ( l - 6 0 - e6f (i + U(1)) 2 , 


(3.19) 


and at infinity of the phase space up) = 0 reduces to the following condition 

w 2 {1) > e(l - 60 + £601 + v {1) ) 2 . (3.20) 

The regions of the phase space where this condition is not fulfilled are nonphysical. 
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Figure 1. The phase space portraits for the dynamical system (3.17) at the invariant 
manifold at infinity defined by = 0. The following values of the model parameters were 
assumed: e = —1, £ < 0 (top left), e = +1, £ < 0 (top right), e = —1, 0 < £ < g (bottom 
left), s = +1, 0 < £ < | (bottom right). 


In figures 1, 2, 3 we presented the phase space diagrams of the system (3.17) at 
the invariant manifold at infinity. The diagrams represent all representative cases for 
all possible values of the non-minimal coupling constant £. The shaded regions denote 
nonphysical regions of the phase space. The dashed lines denote the singularity of the 
time transformation (3.18) where w 2 ^ = e6£(l — On figures 2 and 3 these lines 

are separatrices of the saddle type critical points D\ , D 2 , D 3 , ZD 4 and they intersect 
at the critical point E. 

In table 4 we gathered the eigenvalues of the linearisation matrix for the critical 
points which establish stability conditions for a given asymptotic state. They can be 
used to find the linearised solution in the vicinity of the critical points. The effective 
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Figure 2. The phase space portraits for the dynamical system (3.17) at the invariant 
manifold at infinity defined by up) = 0. The following values of the model parameters were 
assumed: e = —1, g < £ < ^ (top left), e = +1, g < £ < ^ (top right), e = —1, ^ \ 

(bottom left), e = +1, (bottom right). 


equation of state parameters give the physical interpretation of the solutions. 

The second chart of the projective variables Up) leads to the following dynamical 
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Figure 3. The phase space portraits for the dynamical system (3.17) at the invariant 
manifold at infinity defined by = 0. The following values of the model parameters were 
assumed: e = —1, \ < £ < ^ (top left), e = +1, \ < £ < ^ (top right) e = —1, ? > 
(bottom left), e = +1, £ > 2b (bottom right). Note that the presented phase space diagrams 
for the given parameters are topologically equivalent both for canonical and phantom scalar 
fields. The difference comes from dynamical behaviour in the 3-dimensional phase space. 

U(2)(1 + U(2))(W(2) -e6f(l-6f))- 

2 (^(2) + 6£) ^Qa ) o'U( 2 ) + — £ (1 — 6£)it( 2 ) — £ 2£(1 + «(2)) 2 ) , 

u ( 2 )(( m ( 2 ) - 2 )( w ( 2 ) -e6^(l-60) + 

+ 2 (^A,0^ 2 2) + ^ W h) — _ 60^(2) — ^2^(1 + W(2)) 2 j j , 

U(2)W(2)(wf 2) - e6f(l - 60) , 


system 

dtt( 2 ) 

dV(2) 

dV(2) 

dV(2) 

dw( 2 ) 

dr/ (2 ) 


(3.21) 



























“(2) 

w h ) 

existence 

A 

B 

C 

D 

R 

2« 

1-4J 

-6£± V-6£(l-6£) 

-6£ 

0 

0 

0 

±^(1 - 60 
±Ve6£(l ~ 60 

£ < 0 or £ > 7 

e = —1: £<0V£>|or£ = +1: 0 < £ < 7 

e = -l: £> I 

e = —1: £ < 0 or e = +1: £ > 0 


Table 2. Critical points of the system (3.21) in the map 17(2) at infinity of the phase space 
defined by the condition V( 2 ) = 0. 


where the time function along phase space trajectories is 


2 u / 2 

= v, 0 \ — = [w, 


d? 7 ( 2 ) (2) dr 


( 2 ) 


(1 - 60 ) 


din a 


(3.22) 


The critical points of the system at the invariant manifold at infinity given by the 
condition vp) = 0 are gathered in table 2. 

In the projective coordinates of the chart 1/(2) the energy conservation condition 
(2.8) is now 



+ Q 


m, 0 



u ( 2 ) = w p) - - 6 O m ( 2 ) - + «( 2)) 2 , 


(3.23) 


and at infinity of the phase space vp) = 0 it reduces to 

Wp) > e(l - 60w( 2) + £6^(1 + tt( 2 )) 2 • (3.24) 

In figures 4 and 5 we presented the phase space diagrams of the system (3.21) 
at the invariant manifold at infinity. The diagrams represent all representative cases 
for all possible values of the non-minimal coupling constant 0 The shaded regions 
denote nonphysical regions of the phase space where the energy conservation condition 
(3.24) is not fulfilled. The dashed horizontal lines denote the singularity of the time 
transformation (3.22) where iujL = £601 — 6£). 

In table 4 we gathered the eigenvalues of the linearisation matrix of the system 
(3.21) calculated at the critical points. The stability conditions can be directly es¬ 
tablished for different values of the non-minimal coupling constant 0 The values of 
the effective equation of state parameter give the physical interpretation of the given 
asymptotic state. 

The information gathered in table 4 can be used in order to obtain linearised 
solutions in the vicinity of the critical points. In the next section we will concentrate 
our attention on the physical properties of the states. Here we give the linearised 
solutions in the vicinity of the critical points B , which as it turns out, represent very 
interesting asymptotic states from the physical point of view. 

The critical points B in the chart Up) are located at 

(u* {2) = -6£ ± v/-6£(l - 6£), u* 2) = 0 , t^ 2) = o) 
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Figure 4. The phase space portraits for the dynamical system (3.21) at the invariant 
manifold at infinity defined by v^) = 0. The following values of the model parameters were 
assumed: e = — 1 , £ < 0 (top left), e = + 1 , £ < 0 (top right), e = — 1 , 0 < £ < g (bottom 
left), e = + 1 , 0 < £ < | (bottom right). 


and they exist only for £ < 0 or £ > |. The eigenvalues of the linearisation matrix are 

A 1 = -e6£(l~60(3 + 2^ 2) ), 

A 2 = -e6£(l-60(3 + 4 } ), (3.25) 

A 3 = e6f(l - 6 £)w( 2 ) • 

Matching the existence conditions we obtain direct conditions for the stability of the 
given states. For the critical point denoted as B i on the phase space diagrams 4 and 
5 with u * 2 ) = —6£ — y/— 6£(1 — 6£) we have the following conditions for a stable node 
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Figure 5. The phase space portraits for the dynamical system (3.21) at the invariant 
manifold at infinity defined by v^) = 0. The following values of the model parameters were 


assumed: e = — 1 , ^ < f; < (top left), e = +1, i < £ < (top right), e = —1, < £ < 


(middle left), e = +1, 
(bottom right). 


3_ 

16 

<£< 


(middle right) 
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16 

= -!,<£> 


(bottom left), e = +1, £ > 

















































with Xi < 0 


e = —l:£<OV-<£<—, 

0 16 

and for an unstable node with Aj > 0 

e = +l:f<0V^<f<4 . 

6 lo 

The state represented by the critical point B\ is a transient one in the form of a saddle 
otherwise. The critical point denoted as B 2 on the phase space diagrams 4 and 5 with 
u*, 2 , = — 6 £ + yj— 6^(1 — 60 is a stable node with A, < 0 for 


and an unstable node with A* > 0 for 

£ = +1 : £ > x > 
b 

otherwise it represents a transient state. 

Finally, the solutions in the vicinity of those points are in the plain form 

«(2)(b(2)) = u*2) + An (2) e Al?? ( 2 > , 

V(2)(ji{2)) = Ar(2)e A2,(2) , (3.26) 

W(2){r](2)) = Aw( 2 )e A3,,(2) . 


The stability conditions are given for the time ij^) which is directly connected with the 
scale factor via (3.22). In the next section when investigating the physical interpreta¬ 
tion of the given states we will establish stability condition according to the expansion 
of the universe. 

In order to complete dynamical analysis at infinity of the phase space we need to 
investigate the last char of the projective coordinates. The dynamical variables of the 
chart [/( 3 ) lead to the following dynamical system 

dw( 3 ) 


dh(3) 


dV( 3 ) 


= -«(3)(1 

3 


£6^(1-604))- 

|(«(3) + 60(3)) (^A,0W( 3) + \ ~ £ (1 “ 6 0“(3) “ £ 2 £( M (3) + U (3))' 

= n (3) (l-£6^(l-604))> 


0(3) 

dw (3) 3 fry 2 

— = -- Wm (a KoWm 


(3.27) 


1 - e(l - 604 ) - £ 2 0«(3) + V(3)) 2 + 


+ £801-604) 


where the time function along the phase curves is now given by 


= 4 )^ = (1 - £601 - 604 )) 


0(3) 


din a 


(3.28) 
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u (3) 

* 

V (3) 

existence 

M 

0 

0 

V£ g R 

R 

0 


e = —1: £ < 0 or e = +1: £ > 0 

C 

- 6 ^(*3) 

+ 1 

£ = —1: £<0V£>|or£ = +1: 0 < £ < g 

^65(1-65) 

D 


-L 1 

e = - 1: £ > 1 

e(l-4?) 

v^eqi-ej) 


Table 3. Critical points of the system (3.27) in the map U ( 3 ) at infinity of the phase space 
defined by the condition W( 3 \ = 0 . 


In table 3 we gathered the critical points of the system (3.27) located at the 
invariant manifold at infinity, together with conditions for their existence. 

The physical regions of the phase space are determined from the energy conser¬ 
vation condition (2.8). In the projective coordinates of the chart U ( 3 ) we have the 
following 


fo A ,o + ^m,o^j J W( 3 ) = 1 - e(i - 6f)u(3) - e6f(u( 3 ) +%)) 2 (3.29) 

and at infinity of the phase space defined by the condition e 0 we have 

e(l - 60 «( 3 ) + (u( 3) + v {3) ) 2 < 1. (3.30) 

I 11 figures 6 and 7 we presented the phase space diagrams of the system (3.27) 
at the invariant manifold at infinity. The shaded fields correspond to nonphysical 
regions of the phase space. The diagrams represent all the representative cases for 
all possible values of the non-minimal coupling constant £. The dashed lines again 
represent singularities of the time transformation (3.28) where £6£(1 — 6 £)w^ = 1. 
On figure 7 for the phantom scalar held with e = —1 singular lines form horizontal 
separatrices of the saddle type critical points D t . 

In table 6 we give eigenvalues of the linearisation matrix calculated at the given 
critical point together with the effective equation of state parameters. The eigenvalues 
give the stability conditions while the latter informs about physical behaviour at the 
asymptotic state. It what follows we present the linearised solutions in the vicinity of 
the critical points which represent the most interesting states from the physical point 
of view. 

The critical point M with the effective equation of state w e s = 0 located at 
(«(3) = 0, < 3 ) = 0, «fo = 0), 

exists for any value of the non-minimal coupling constant £ and has the following 
eigenvalues 

Ai, 2 = —j(3±y3(3-16£)), A s = |. 
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Figure 6. The phase space portraits for the dynamical system (3.27) at the invariant 
manifold at infinity defined by uy 3 ) = 0. The following values of the model parameters were 
assumed: e = — 1, £ < 0 (top left) and e = +1, £ < 0 (top right), e = — 1, 0 < £ < | (bottom 
left) and e = +1, 0 < £ < g (bottom right). 

In any possible case it represents a saddle type critical point. The linearised solutions 
are 

U ( 3 )(V( 3 )) = v / 3 ( 3 2 - 16 ^) (( 3 ^ At;(3) “ A i AM (3)) eAlT,(3) - ( 3 £Au (3) - A 2 A u ( 3 ) )e X2 ^ , 
V ( 3 )M = ~ AU(3) ^ Al,?(3) “ ( AlAu ( 3 ) ~~ Ait (3) )e A2??(3) ^ , 

3 

W( 3 )(V( 3 )) = A w { 3 ) e^) , 

(3.31) 

where Art( 3 ) = ufy, Au( 3 ) = and Au>( 3 ) = are the initial conditions. 
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Figure 7. The phase space portraits for the dynamical system (3.27) at the invariant 
manifold at infinity defined by wr 3 ) = 0. The following values of the model parameters were 
assumed: e = — 1 , | (top left), e = + 1 , £ < £ < (top right), e = — 1 , £ > ^ 

(bottom right), e = + 1 , £ > ^ (bottom right). 


The critical points R with the effective equation of state parameter w e s = | are 
located at 


«(3) = °> v k = ± \l^> w h) = 0 


They exist only for e£ > 0, i.e., for the canonical scalar field with e = +1 only for 
£ > 0 while for the phantom scalar field with e — — 1 only for £ < 0. The eigenvalues 
are 


Ai = -6£, A 2 = 6£, A 3 = 12£, 


and for any possible values of the non-minimal coupling constant £ it represents a 


23 



























saddle type critical point. The linearised solutions are 


«(3)fa(3)) = \ (Au (3) - Au (3) ) e 6 ^ (3) + i (Am (3) + Au ( 3) ) e 6 ^ 3 > , 

W(3)(h(3)) = v* {3) - i (Au (3) - Au ( 3} ) e _6 ^ (3) + i (Au (3 ) + Au (3) ) , ^ 3 ’ 32 ) 

W(3)(V(3)) = Aw ( 3) e 12 ^ (3) , 


where Arq 3 ) = Av@) = — v*.^ and Au>( 3 ) = are the initial conditions. 

Note that this critical point denoted on the all relevant phase space diagrams 
as R in the original variables corresponds to the condition e^n 2 (j) 2 = 1, which in the 
Einstein frame formulation of the theory leads to vanishing of the conformal factor and 
the conformal metric degenerates. For an isotropic evolution this is a regular point 
while a strong anisotropic curvature singularity forms there as it was show in [79] (see 
also [80]). In the Jordan frame, followed in this paper, this condition is represented 
by the critical point of a saddle type where two of the separatrices form the boundary 
between physical and nonphysical regions of the phase space. 

The last two critical points which are interesting form the physical point of view 
are denoted on the phase space diagrams as C\ and C 2 . The effective equation of state 
parameter calculated at these points give an infinite value pointing to the presence of 
singularity. The coordinates of the asymptotic states under considerations are 



-6^(3) , V( 3 ) = ± 


V^6£(l - 60 ’ 


W *3 ) = 0 


with the eigenvalues 

A 1 = A 3 = 6£, A 2 = 12£ . 

Together with the existence condition from table 3 we have that the critical point is in 
the form of a stable node only for the phantom scalar held with £ — —l and for £ < 0 
and is an unstable node for the phantom scalar held for £ > | while for the canonical 
scalar held £ = +1 it exists and is in the form of an unstable node for 0 < £ < The 
linearised solutions are 


U(3)(V(3)) = U(3) + (Au ( 3) + 2(1 - 30Au (3) ) e 6 ^ 3 > - 2(1 - 3£)Au (3) e 12 ^ 3 > , 

V(3)(V(3)) = U(3) + A 17 ( 3 ) e 12 ^ 3 ) , (3.33) 

W(3)(V(3)) = Aw ( 3 ) e 6 ^ 3) , 


where Atq 3 ) = u^l — v *.^, Aw( 3 ) = vQ ~ v * 3 ) and Aiu( 3 ) = are the initial conditions. 

What remains is a discussion the other non-hyperbolic critical points with one 
eigenvalue vanishing. In order to obtain a qualitative behaviour in the vicinity of such 
critical points one is forced to use the centre manifold theorem [22, 74, 75]. This is 
beyond the scope of the current study and will be presented elsewhere. 
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Table 4. The eigenvalues and the effective equation of state parameter for the critical points 
at infinity of the system (3.17) in the map gathered in table 1. 



eigenvalues 

W e s 

A 

B 

C 

D 

R 

( £ 3g (1 - 6 f_ ( ^ 160 , -e95(l -6f), -£l25 2 ±5§§) 

( - £65(1 - 60(3 + 2 U* 2) ), -£65(1 - 65)(3 + u* 2) ), £65(1 - 65)^ 2) )) 
(e365 2 (1 - 65), £365 2 (1 - 65), e725 2 (1 - 65)) 

(£125(1 -65)1^, 0, -e125(1 -65)^ 2) ) 

( - 65v/i65™ ( * 2) , e 725 2 , 65v/i65™ ( * 2) ) 

4S 

3(1—4?) 

i+H) 

±oo 

1 

3 


Table 5. The eigenvalues and the effective equation of state parameter for the critical points 
at infinity of the system (3.21) in the map [/( 2 ) gathered in table 2. 

4 Physics from dynamics 

In the previous section we obtained the linearised solutions of the dynamics in the 
vicinity of the stationary states, i.e. critical points. Now we want to extract a physical 
meaning from these solutions in order to find a physical interpretation of specific states 
of the universe represented by the critical points. 

In tables 4, 5 and 6 we have gathered the effective equation of state parameters and 
the eigenvalues of the linearisation matrix for all the critical points in all three charts 
covering the phase space at infinity. The eigenvalues represent stability conditions for 
a given state in the generic case, i.e. for a generic value of the non-minimal coupling 
constant £. 

We start our discussion from the two critical points in the chart 17(3) which rep¬ 
resent a barotropic matter domination and a radiation dominated epoch in evolution 
of the universe. 

The critical point with the effective equation of state parameter w e g = 0 denoted 
on the phase space diagrams as M has at least one unstable direction. From the 
coordinates transformation in the chart f/( 3 ) (3.16) we have that 

H(a) _ 1 

H(a 0 ) w (3 )(77 (3) ) ’ 

and using the last linearised solution form (3.31) we have that the Hubble function in 
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the vicinity of this critical point behaves as 


H(a) 1 3 

—— =- e _ 2 ^ 3 > . 

H(a 0 ) A w( 3 ) 

The time transformation (3.28) up to linear terms in initial conditions gives 

K ln (J>) ■ 

where a is the scale factor and is the initial value of the scale factor in the vicinity 
of the critical point. Finally we obtain that the Hubble function as the function of the 
scale factor is 

~ ff(oW) /_a_ r ! 

H(a o) H(a 0 ) W-V ’ { ' 1 

which is exactly the Hubble function in the barotropic dust matter domination phase 
of expansion of the universe. 

The asymptotic states denoted on the phase space diagrams as R have the effective 
equation of state parameter = |. For arbitrary values of the non-minimal coupling 
constant £ they are in the form of a saddle type critical points. Again from the 
coordinates of the chart (3.16) we have that 

H(a) _ 1 

H(ao) w ( 3 )( 77 (3) ) ' 


From the linearised solutions in the vicinity of this state (3.32) we have the Hubble 
function 


H(a) 


e - 12 b7 ( 3) 


H(a 0 ) Aw (3) 

The time transformation (3.28) up to linear terms in initial conditions leads to 


din a ~ ( 6 £ + e 6 f(l - 6 f)u* 3 )((Au (3 ) - Au (3) )e 6?T?(3) - (Au (3) - Av^e 6 ^) )d 7 ? (3) , 


integrating this expression and expanding for small values of the time r/( 3 ) we obtain 


111 (^)) ~ 6 ^ 3 )( 1 “ £2 ^ ~ 6 0 ^( 3 ) A ^( 3 )) , 

where, again, a is the scale factor and a® is the initial value of the scale factor in the 
vicinity of the critical point. The resulting Hubble function is 

H(a) H W / a \~ 2 f a \ 

W{ao) ~ wo J wo J ’ { } 

which is the Hubble function during the radiation domination epoch in the history of 
the universe with the small correction proportional to initial conditions and the value 
of the non-minimal coupling constant. 
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eigenvalues 

Weft 

M 

( _ 1 (3 + ^3(3 - 16?)) , (3 - ^3(3 - 16$)) , f) 

0 

R 

( - 6£, 65, 125) 

1 

3 

C 

(65, 125, 65) 

zb OO 

D 

(e!25(l - 65)^3)^), -e 12 ^ 1 “ e i) u (3) v h)’ °) 

- 


Table 6. The eigenvalues and the effective equation of state parameter for the critical points 
at infinity of the system (3.27) in the map [/( 3 ) gathered in table 3. 


Investigating dynamics of cosmological models within various theories of gravity 
one usually encounters some sort of singularities. They can be in the form of an initial 
finite scale factor singularity [87, 88 ], future finite scale factor singularities: the sudden 
future singularities [89], the Big Brake singularity [90] or the Big Boost singularity 
[91]. The next two asymptotic state represent precisely a finite scale factor singularity. 

The critical points denoted on the phase space diagrams as C\ and C 2 in the chart 
f/( 3 ) have the following coordinates at the invariant manifold at infinity W ( 3 ) = 0 



-6^(3) , V 


(3) 


Veeea-eo 


are in the form of a stable or a unstable node depending on the sign of the non-minimal 
coupling constant £. Using the linearised solution (3.33) we obtain the following form 
of the Hubble function 


H(a) _ 1 

H(ao) W(3)(ri(3)) 


_ e -6b?( 3 ) _ 

Aw (3 ) 


Clearly for £ > 0, i.e. when the asymptotic state is in the form of an unstable node, 
when approaching this state ??( 3 ) —* — 00 the Hubble function goes to infinity. The 
same is true for £ < 0 , i.e. when the state is in the form of a stable node and at 
the limit 77 ( 3 ) —>■ + 00 . Now we want to show that the scale factor at this state has a 
finite value. From the definition of time 77 ( 3 ) (3.28) and up to linear terms in initial 
conditions we have the following 


din a « —e: 12 £(l - 6 £)u ( * 3) Au ( 3 ) e 12 ^ 3 >d 77 (3) , 

and after integration we obtain 

/ a (/) \ 

ln ~ £<<1 ~ 6 ^ v h A m > 


where a W and are the initial and the final values of the scale factor. Interestingly, 
the finite scale factor singularity under considerations can appear in two types. Namely 
in one case when > a® and in opposite case when < aS l \ it means that it can 
be either in the future as well as in the past according to the expansion of the universe. 
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Which scenario is realised crucially depends on the form of the scalar field as well as 
the value of the non-minimal coupling constant £. For example in the bottom right 
phase space diagram on figure 6 asymptotic state represented by the point C\ for the 
canonical scalar field £ — +1, 0 < £ < | and initial conditions taken in the physical 
domain (white region) we have e(l — 6^)v* 3 ^Avp) < 0 and hence < aP\ the finite 
scale factor singularity lies in the past according to the expansion of the universe. The 
same reasoning can be applied to the phantom scalar field £ — — 1 in the top left phase 
space diagram on figure 6 and the phase space diagrams in the left panel of figure 7. 

The asymptotic states in the form of critical points discussed so far were previously 
found for a wide class of the scalar field potential functions [22], Now we can move on 
to investigations of the physical properties of the critical points previously missed due 
to lack of the dynamical analysis at infinity of the phase space. 

Two critical points in the charts Up) and Up) denoted on the phase space diagrams 
as B i and B 2 possess very interesting and special properties. Namely for the specific 
value of the non-minimal coupling constant £ = A the effective equation of state 
parameter calculated at those points are Wes\ Bi = — 1 and w e ,g| = 0. We start from 

the latter one. The linearised solutions (3.26) specified for the critical point B 2 and 
for £ = A; are 

3 27 

U( 2 ) 4 ( 2 )) = _ 4 + , 

V(2)(V(2)) = Av( 2 )e £ i^(2) ; ( 4 -3) 

27 

W(2){V(2)) = A wp)e £ ^ v ^ , 


where we have a stable node for the phantom scalar field e — — 1 and an unstable node 
for the canonical scalar field e = +1. 

From the coordinates of the chart Up) (3.16) and using those linearised solutions 
we have 

H ( a ) = = H) e - £ §, (2) 

H(a 0 ) vp)(-np)) Au (2) 

and from (3.22) up to linear terms in initial conditions 

in (Jr) k 4’« - 

Finally the Hubble function in the vicinity of the asymptotic state is given by 

H(a 0 ) H(a 0 ) 

It represents unstable during expansion the Einstein-de Sitter solution for both types 
of the scalar field. 

The critical point B\ with the effective equation of state parameter w e s\ B = — 1 
for £ = 4 has the vanishing eigenvalue in the direction of the variable up) and is of a 
non-hyperbolic type. Its detailed discussion deserves a special treatment and will be 
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presented in the next section. Nevertheless, based on the linearised solutions in the 
remaining directions we can find the behaviour of the Hubble function in the vicinity 
of this point. From the last two linearised solutions of (3.26) and for £ = A we have 
the following 

v ( 2)M = Av (2) e e ^m , 

W(2)(V(2)) = A W(2)e £ ^ V{2) , 


and from the coordinates of the chart U ( 2 ) (3.16) we directly find that the Hubble 
function 

H ( a ) _ W (2)(^(2)) = AU( 2) = 

H(a 0 ) v {2) (r, {2) ) A w {2) 

is constant during the evolution in the vicinity of this state and corresponds to the de 
Sitter type of evolution. Because of the non-hyperbolicity of this state we are unable 
to determine stability conditions on the basis of the linearisation procedure and we are 
forced to take a different path in order to find specific terms leading to the unstable 
de Sitter state during expansion of the universe. 


5 Unstable de Sitter state at the bifurcation value £ = ^ 


The asymptotic states of dynamical systems of cosmological origin correspond to dif¬ 
ferent phases of evolution of universe. In analysis presented so far we found that there 
is a special bifurcation value of the parameter £ for which the de Sitter state exists 
at infinity of the phase space. The critical point corresponding to this state is of a 
non-hyperbolic type with one vanishing eigenvalue. I 11 order to find solutions in the 
vicinity of critical points of such type one can use the standard approach, the centre 
manifold theorem [22, 74, 75]. We will take a different approach to emphasise and 
extract the physical significance of resulting dynamical behaviour in the vicinity of 
this critical point. 

We begin with the system (3.21) together with the following transformation of 
the phase space variables 


« = '“(2) + 7. , V = “(2) . W = “(2) . 


(5.1) 


and the dynamical system is 


d u 

<ty(2) 

dh 

<ty(2) 

dw 

d??(2) 


= - \u - 


= — 2v \ \ u — 


u — 


q \ q 

™ + £ m)-2' U - 


9 


„ _ 1 1 2 3 

^A,ov + -w - e-u + e— 
3 4 16 


^ + £ 64 ) + b( nA .°^+ + e 


16 


= —2 wii 


a --ll® + e| E 


(5.2) 
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where 


d _d /_ 9 

drj (2) " "d7 ~~ V l,+£ ei4/dlna 


(5.3) 


Now the critical point under investigations is located at the centre of the coordinate 
system and we were able to lower the degree of dynamical equation in variables v and 

w. 

Up to linear terms in the variables we obtain the following dynamical system 


9 9 _ 

U « — ilAflV ~ —W, 
16 16 

27 _ 

27 

W ~ £ - W , 

64 ’ 


(5.4) 


where the prime denotes differentiation with respect to the 77 ( 2 ) time and the last two 
equations are now decoupled. This linearised system can be easily solved, resulting in 


u(V( 2 )) = u (l) - - W w ) ^1 - exp ^£^77(2)^ , 


v(V( 2 )) = " W 


27 


exp l £ ^( 2 ) I ’ 


(5.5) 


Mm) = exp [ £—?7( 2 ) 


27 


We need to note that the critical point corresponding to the de Sitter state is located 
at the centre of the coordinate system. Above solutions for a generic set of initial 
conditions do not lead to this point. We need to take some open and dense subset of 
possible initial conditions satisfying 


u 


w = e 5 (!2A,o® (i) - i5 (i) ) • 


(5.6) 


The effective equation of state parameter up to linear terms in the phase space variables 

is 

1 ^AfiV + _ i ^A,0^ — 


^eff 






(5.7) 


and inserting above solutions and taking the limit £77(2) —> — oo we obtain pure de 
Sitter state. From the other hand the energy conservation condition 


Q 


m, 0 


a 0 


-3 


1 9 

= h 2 - U A ,o + ~ s-{x + z) 2 > 0 


(5.8) 


calculated for the specific initial conditions (5.6) indicate that the effects of the baro- 
tropic dust matter must be negligible in order to obtain the de Sitter state in form of an 
unstable node. In the upcoming paper [92] we investigate in details stability conditions 
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and asymptotic forms of scalar field potential functions giving rise to a singular free 
type of cosmological evolution. On figure 8 we presented the phase space portrait of 
the dynamics of the model for the canonical scalar field £ — +1 without the barotropic 
matter and £ = There is open and dense subspace of initial conditions leading to 
a non-singular evolutional path from the unstable de Sitter state to the stable one. 

The obtained possible evolutional path can indicate for some generalisations of 
the seminal Starobinsky type of evolution [93], as well as obtained evolution was not 
expected nor designated from the very beginning [94, 95] and may suggest some addi¬ 
tional symmetry in the matter sector of the theory. 



Figure 8. The phase space portrait compactified with a circle at infinity of the model 
without the barotropic matter and the non-minimal coupling constant £ = The de Sitter 
state at the centre dS_|_ is a global attractor for the dynamics. Four critical points at infinity 
describe unstable de Sitter states dS_ and Einstein-de Sitter states EdS_. All the trajectories 
connecting two de Sitter states correspond to possible singularity free evolutional paths. 


6 Beyond linear approximation 

Dynamical systems theory gives us great opportunity to investigate qualitative be¬ 
haviour of dynamical systems in vicinity of stationary states. The linearisation proce¬ 
dure applied to hyperbolic critical points enables us to obtain description of dynamics 
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up to some characteristic time interval assuming some initial conditions. We were us¬ 
ing this approach extensively in the previous sections. Unfortunately sometimes it is 
not enough to recover physical description of the processes under investigations. We 
need to go beyond the linear approximation. 

Within local dynamical systems theory in addition to the linearisation procedure 
we are equipped with the notion of invariant manifolds [74, 75]. 

In section 2 we stated that the system (2.11) is equipped with two invariant man¬ 
ifolds which reduce dimensionality of the system. One invariant manifold corresponds 
to the static Einstein universe with the vanishing Hubble function and the second one 
corresponds to the de Sitter universe with the constant value of the Hubble function. 
In what follows we use the latter in order to obtain the exact solutions of reduced 
dynamics and then we obtain the Hubble function beyond the linear approximation. 

Let us begin with the dynamical system (2.11) in full form 



( 6 . 1 ) 


din a 


H 2 ’ 

H 


where now we treat h 2 as the phase space variable and the acceleration equation is 


H _ o 3 11 a , o + \fi 2 - e(l - 6£)x 2 - e2£(x + z) c ‘ 
lp~~ + 2 ft 2 - e6£(l - 60s 2 


( 6 . 2 ) 


On the invariant manifold = 0 corresponding to the de Sitter state the dy¬ 
namical system (6.1) reduces to a simple linear system 


dx 

din a 
dz 

din a 


= — 3x — 12 £z , 


= x. 


(6.3) 


This linear system of equations can be directly integrated. We make the following 
change of variables 


x = Ai p + A 2 q , 
z = p + q , 


(6.4) 


where Ai — —|— |^/l - A 2 — 
matrix of the linear system (6.3). 


—| + § y 1 — are the eigenvalues of the constant 
Then we obtain 


dp 

din a 
d q 

din a 


AiP, 


A 2 <?- 


(6.5) 
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and the solutions are 


p(a) = p(a 0 ) 


q(a) = q(a 0 ) 


a o 






^2 


( 6 . 6 ) 


These solutions are the exact solutions of the reduced system, thus the initial conditions 

can be arbitrary chosen and we choose them at the present epoch. The quantities p(a 0 ) 

and q(a 0 ) denote values of the phase space variables at a 0 — 1 value of the scale factor. 

The third equation in the system (6.1) expanded in to the Taylor series up to 

linear terms in z variable is 

d h 2 / \ 

— « 3 fi A , 0 - 3 h 2 - £3 ((1 - 4£)x 2 + 4£xz) , (6.7) 

where dN = din a. 

The condition z 2 ~ 0 expressed in the variables p and q leads to 


z 2 = (p + q) 2 = p 2 + 2 pq + q 2 


0. 


For the remaining variables appearing in equation (6.7) we obtain 

2 


9 


x = t (p + q) + \ 1 - (p ~ q) 


16 




( 6 . 8 ) 


and 


xz = 


(p + q) + 


16 


€(p-q) (p + q ) 


16 


f ( p 2 - q 2 ) ■ (6-9) 


Finally, inserting those approximation in to the equation (6.7), the dynamical equation 
for the Hubble function is given by 

« —3 h 2 + 3f2 Ai0 + ^ 9 ^1 —— ^(Ai + 6£)p(lV) 2 + (A 2 + 6£) g(fV) 2 j , (6.10) 

where the functions p(N) and q(N) are the exact solutions of the system (6.3) reduced 
to the invariant manifold of the de Sitter state and are given by (6.6). 

Now we are ready to integrate the dynamical equation (6.10) and the result is 


(S) a “ n « +n 


m, 0 


ao 


-3 




-3 


(Ai + 6£) p^ 


Oq 
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( 6 . 11 ) 


+ (A 2 + 6£) q 0 ( — ] 
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The first two terms in this solution correspond to the standard ACDM model and 
remaining term can be interpreted as an extension depending on the non-minimal 
coupling £ and the current values of the phase space variables p(ao) and q(ao). For 
£> ^ this additional term reduces to an oscillatory behaviour with dumping. We need 
to note that in the first approximation beyond linear investigated in this section, the 
obtained Hubble function reduces to the ACDM model with the assumption £ = 

7 Conclusions 

In this paper the global dynamics of the simplest case of the non- minim a,lly coupled 
scalar field cosmology was investigated in the details in the framework of dynamical 
system methods. The dynamics is simple in this sense that in the flat Friedmann- 
Robertson-Walker model with the dust barotropic matter can be reduced to the three- 
dimensional dynamical system. The complexity of the dynamics is revealed while 
considering different values of the non-minimal coupling constant £. The bifurcation 
values of this parameter were identified and the global phase analysis for representative 
cases was performed. 

It is interesting that dynamics with the minimal £ = 0 and the conformal £ = | 
coupling constants are equivalent to the ACDM model with additional terms which 
behave like a stiff matter and a radiation matter, respectively. In the case of conformal 
coupling for the canonical scalar field £ — +1 this additional term behaves like an 
ordinary radiation while for the phantom scalar field £ — —1 like a dark radiation with 
a negative energy density. The corresponding Hubble functions were obtained in the 
exact form are the models are therefore integrable by elliptic functions. 

We demonstrated that for the case of the vanishing potential function of the scalar 
field the right hand sides of the dynamical system are homogeneous functions of the 
third degree which gives the possibility of lowering the dimension of the phase space 
by one using the projective coordinates. 

We have shown that the constant £ describing the non-minimal coupling between 
the curvature and the scalar field is the crucial ingredient not only during the inflation¬ 
ary epoch [17] but also is describing the current accelerated expansion of the universe. 
While the special bifurcation value of the non-minimal coupling £ = 3/16 is distin¬ 
guished by the dynamics of the model and can lead to a singularity free cosmological 
evolution may suggest some additional symmetry in the matter sector of the theory. 
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